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Use  of  X-Ray  Diffraction,  Molecular  Simulations,  and 
Spectroscopy  to  Determine  the  Molecular  Packing  in  a 
Polymer-Fullerene  Bimolecular  Crystal 


Nichole  Cates  Miller,  Eunkyung  Cho,  Matthias  J.  N.  Junk,  Roman  Cysel,  Chad  Risko, 
Dongwook  Kim,  Sean  Sweetnam,  Chad  E.  Miller,  Lee  J.  Richter,  R.  Joseph  Kline, 
Martin  Heeney,  lain  McCulloch,  Aram  Amassian,  Daniel  Acevedo-Feliz, 

Christopher  Knox,  Michael  Ryan  Hansen,  Dmytro  Dudenko,  Bradley  F.  Chmelka, 
Michael  F.  Toney,  Jean-Luc  Bredas,  and  Michael  D.  McGehee'k 


Bullc-heterojunction  (BHJ)  solar  cells  made  from  blends  of 
semiconducting  polymers  and  fullerenes  are  attractive  because 
they  can  be  printed  at  low  cost  and  have  demonstrated  effi¬ 
ciencies  greater  than  8%.I1-5!  For  these  solar  cells  to  produce 
current,  excitons  created  during  light  absorption  must  reach 
a  conjugated  polymer  (donor):  fullerene  (acceptor)  interface 
and  dissociate  by  charge  transfer.  Electrons  and  holes  must 
then  travel  through  the  fullerenes  and  polymer,  respectively, 
to  reach  the  electrodes  before  recombination  occurs.  The  for¬ 
ward  and  back  charge-transfer  processes  and  charge  transport 
depend  critically  on  the  molecular  packing.  Slightly  modifying 
the  distance  between  the  donor  and  acceptor  can  change  the 
charge-transfer  rates  by  more  than  an  order  of  magnitude  and 
consequently  have  an  enormous  impact  on  exciton  dissociation 
and  recombination.!6!  Likewise,  small  changes  in  the  molecular 
packing  can  significantly  change  charge-carrier  mobilities. !7,81 
Although  important  processes,  such  as  charge  separation  and 


recombination,  can  depend  on  the  polymer-fullerene  wave- 
function  overlap,  very  little  is  known  about  how  polymers  and 
fullerenes  pack  on  the  molecular  level  since  it  is  very  difficult  to 
accurately  characterize  molecular  packing  in  BHJs. I9' 

To  design  materials  for  BHJs  rationally,  it  is  therefore  critical 
to  develop  an  understanding  of  the  molecular  packing  and  its 
effects  on  key  electronic  processes.  Until  now,  it  has  not  been 
possible  to  establish  the  detailed  molecular  organization  of 
BHJ  blends  using  X-ray  diffraction  (XRD),  because  the  blends 
were  either  amorphous  or  too  disordered  to  yield  enough  XRD 
reflections  to  accurately  determine  the  unit  cell  and  molecular 
packing.  The  best  BHJ  material  combinations  have  largely 
been  discovered  by  trial  and  error.  Currently  no  one  is  able  to 
accurately  predict  why  certain  material  combinations  perform 
dramatically  better  than  others  or  why  certain  processing  con¬ 
ditions  can  improve  performance  in  some  materials  combina¬ 
tions,  yet  reduce  performance  in  others. 
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A)  B) 


Figure  1.  The  chemical  structures  of  PBTTT  (A)  and  PC71BM  (B).  R  can 
be  C12H25,  C14H29,  or  C16H33. 


Here,  we  determine  for  the  first  time  the  detailed  struc¬ 
ture  of  a  BHJ  blend  that  establishes  at  the  molecular  level  the 
properties  and  function  of  BHJ  solar  cells.  Specifically,  we 
use  a  combination  of  XRD  techniques,  molecular  mechanics 
(MM)  and  molecular  dynamics  (MD)  simulations,  2D  solid- 
state  NMR  spectroscopy,  and  IR  spectroscopy  to  determine  the 
unit  cell  and  atomic  positions  of  the  poly-(2,5-bis(3-alkylthi- 
ophen-2-yl)thieno[3,2-b]thiophene)  (PBTTT) :phenyl-C71-butyric 
acid  methyl  ester  (PC71BM)  bimolecular  crystal,  in  which  the 
PC71BM  molecules  intercalate  between  the  PBTTT  side  chains 
(see  Figure  1  for  molecular  structures).!10!  We  find  that  the 
interactions  between  PBTTT  and  PC71BM  are  strong  enough 
to  severely  disrupt  the  PBTTT  packing  by  inducing  bends 
and  twists  in  the  PBTTT  backbone,  as  well  as  bends  in  the 
side  chains,  all  of  which  allow  increased  interaction  between 
the  PBTTT  and  the  PC71BM.  This  finding:  i)  has  significant 
implications  for  other  polymeufullerene  blends,  because  strong 
polymer-fullerene  interactions  can  cause  polymer-fullerene 
intermixing  in  both  semi-crystalline  and  amorphous  blends; 
and  ii)  is  surprising  considering  that  until  recently  BHJs  were 
thought  to  contain  phase  separated  regions  of  relatively  pure 
polymer  and  pure  fullerene  with  very  little  mixing  at  the  molec¬ 
ular  scaleJ11"14'  Our  work  lends  insight  into  the  local  packing  in 
other  blends  with  fullerenes  present  in  the  amorphous  regions 
of  semi-crystalline  polymers.  The  methodology  presented  here 
makes  it  possible  to  determine  the  details  of  how  polymers  and 
fullerenes  are  arranged  relative  to  each  other. 

We  recently  generated  a  phase  diagram  for  PBTTT- 
C14:PC71BM  blends!15!  and  demonstrated  that  the  highest- 
efficiency  PBTTT-C14:PC71BM  solar  cells  contain  80  wt% 
PC71BM.!1S'17!  At  this  composition,  approximately  half  of  the 
film  consists  of  PBTTT-C14:PC71BM  bimolecular  crystals  in 
which  long-distance  hole  transport  occurs  and  the  other  half 
consists  of  regions  of  fullerene  that  allow  long-range  electron 
transport.  As  the  PC73BM  fraction  of  the  film  decreases  from 
80  to  50  wt%,  the  power  conversion  efficiency  drops  from 
2.35  to  0.16%  due  to  the  elimination  of  percolating  fullerene 
clusters  that  are  needed  to  provide  electrons  with  adequately 
fast  pathways  to  the  electrodes. 

PBttt:PC71BM  bimolecular  crystals  form  upon  spin¬ 
coating  blends  of  PBTTT  and  PC7jBM  from  a  variety  of  sol¬ 
vents,  as  well  as  upon  thermal  annealing  or  solvent  annealing 
of  PBTTT:PC71BM  bilayer  films  (Supporting  Information, 
Figure  SI).  The  structure  of  these  bimolecular  crystals 
appears  to  be  independent  of  the  solvent  and  processing  con¬ 
ditions  as  long  as  a  brief  anneal  is  provided  to  allow  the  crys¬ 
tals  to  approach  thermodynamic  equilibrium.  Diffraction 
from  both  the  PBTTT-C14:PC71BM  bimolecular  crystal  and 


pure  PBTTT-C14  is  observed  in  blends  with  30  and  40  mol% 
PC71BM,  indicating  that  these  blends  contain  excess  polymer 
(Figure  2  A).  Blends  with  60  mol%  PC71BM  exhibit  diffraction 
from  the  bimolecular  crystal  and  pure  PC71BM,  indicating  that 
these  blends  have  excess  PC71BM.  Diffraction  from  only  the 
bimolecular  crystal  is  observed  in  blends  with  exactly  50  mol% 
PC71BM  (=60  wt%  PC71BM),  which  corresponds  to  one  PC71BM 
molecule  per  PBTTT-C14  monomer.  The  purpose  of  this  study 
is  to  determine  the  detailed  molecular  structure  of  the  PBTTT- 
C14:PC71BM  bimolecular  crystal.  Consequently,  all  of  the  exper¬ 
iments  described  in  the  remainder  of  this  manuscript  were 
performed  on  films  with  50  mol%  PC71BM  so  that  the  pure 
bimolecular  crystal  could  be  studied.  Although  these  films  do 
not  make  efficient  solar  cells,  the  structural  insights  obtained 
are  directly  related  to  key  molecular  interactions  and  processes 
that  govern  macroscopic  photovoltaic  performance. 

Figure  2B  shows  the  experimental  2D  grazing  incidence 
X-ray  scattering  (2D  GIXS)  pattern  for  the  PBTTT-C14:PC71BM 
bimolecular  crystal.  The  intense  (hOO)  reflections  along  the 
nominally  qz  direction  demonstrate  the  formation  of  a  lamellar 
structure  with  the  side  chains  oriented  out  of  the  plane  of  the 
substrate.  As  typically  found  for  spin-cast  films,  the  film  is  iso¬ 
tropic  in  the  substrate  plane,  making  it  difficult  to  distinguish 
between  (OkO),  (00!)  and  (0k!)  reflections,  which  all  occur  near 
the  axis.  Hence,  two  sets  of  unit-cell  parameters  are  compat¬ 
ible  with  the  limited  number  of  experimentally  available  reflec¬ 
tions.  Both  the  orthorhombic  lattice  parameters  a  =  31  A,  b  = 
12.8  A,  and  c  =  13.5  A  and  the  monoclinic  unit-cell  parameters 
a  =  31  A,  b  =  10  A,  c  =  13.5  A,  and  a  =  108°  index  all  of  the 
observed  GIXS  reflections.  To  distinguish  between  (0k!)  reflec¬ 
tions  and  identify  the  correct  unit-cell  parameters,  we  prepared 
a  bimolecular-crystal  film  with  preferred  crystallographic  ori¬ 
entation  parallel  to  the  substrate  by  aligning  pure  PBTTT-C14 
using  a  flow-coating  method,!18!  spin-coating  a  PC71BM  layer 
on  top  of  the  aligned  PBTTT-C14  film  using  the  orthogonal 
solvent  dichloromethane,  and  annealing  the  resulting  film  to 
allow  diffusion  between  the  layers  and  the  formation  of  aligned 
bimolecular  crystals  (Supporting  Information,  Figure  SI).  The 
in-plane  angular  dependence  of  the  reflection  intensities  of 
the  aligned  film  at  q ^  =  0.49  A-1  and  q^  =  0.66  A-1  shows  that 
the  correct  unit  cell  is  monoclinic,  rather  than  orthorhombic 
(Supporting  Information,  Figure  S2,S3). 

Solid-state  NMR  spectroscopy  is  sensitive  to  local  composi¬ 
tional,  structural,  and  dynamical  environments,  in  both  ordered 
and  disordered  regions  of  the  materials.!19,20!  Solid-state  NMR 
spectroscopy  therefore  provides  complementary  insight  on  the 
molecular  structures  of  PBTTTfullerene  blends,!21!  including 
disordered  moieties  for  which  XRD  is  less  useful.  In  particular, 
the  solid-state  2D  13C{1H}  heteronuclear  correlation  (HETCOR) 
NMR  spectroscopy  technique  provides  detailed  information 
about  the  molecular  interactions  between  PBTTT-C16  and 
PC71BM  by  identifying  dipole-dipole-coupled  13C  and  3H  nuclei 
that  are  associated  with  molecularly  proximate  moieties.  For 
example,  Figure  3A  shows  a  2D  13C{1H}  HETCOR  spectrum 
acquired  for  a  PBTTT-C16:PC71BM  blend  with  a  molar  ratio 
of  1:1.  This  spectrum  exhibits  strong  2D  intensity  correlations 
(red,  blue  arrows)  between  13C  and  3H  signals  from  the  conju¬ 
gated  polymer  (orange)  and  the  fullerene  derivative  (blue)  that 
establish  their  mutual  molecular  proximities  (<1  nm).  Such 
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Figure  2.  A)  Specular  XRD  patterns  for  PBTTT-C14:PC71BM  blends  with  a  range  of  molar  ratios.  B)  Experimental  2D  GIXS  pattern  for  the  bimolecular 
crystal.  The  highest-intensity  non-specular  reflections  are  circled  for  comparison  to  the  simulated  pattern  in  (D).  C)  Specular  XRD  reflection  intensi¬ 
ties  normalized  by  the  (TOO)  peak  intensity.  The  experimental  data  are  shown  in  blue,  and  the  intensities  calculated  from  the  structure  factor  of  the 
simulated  unit  cell  are  shown  in  orange.  D)  Calculated  2D  GIXS  pattern  for  the  simulated  unit  cell.  A  list  of  the  diffraction  peak  positions  is  given  in 
Table  S2  in  the  Supporting  Information. 


intensity  correlations  provide  unambiguous  proof  of  the  inter¬ 
calation  of  the  fullerene  derivative  in  the  polymer  network. 

Specifically,  signals  associated  with  13C  moieties  from  the 
PC71BM  side  group  (e.g.,  at  172  ppm  from  the  C=0  groups,  at 
129  ppm  from  the  phenyl  rings,  and  at  36  and  73  ppm  from  the 
cyclopropane  anchor)  are  correlated  with  1H  signals  at  1.4  ppm 
from  the  aliphatic  PBTTT-C16  polymer  side  chains.  Addi¬ 
tional  intensity  correlations  are  observed  between  13C  signals 
from  the  C7o  fullerene  moieties  and  3H  signals  at  1.4,  7.2,  and 
8.7  ppm  from  the  aliphatic  side  chains  and  the  aromatic  back¬ 
bone  protons  of  PBTTT-C16,  respectively.  These  correlated  2D 
NMR  spectroscopy  signal  intensities  establish  strong  interac¬ 
tions  and  hence  close  molecular  proximities  between  these 
molecular  moieties.  Furthermore,  separate  solid-state  1H{1H} 


double-quantum  NMR  spectroscopy  analyses  enabled  the 
local  distances  between  different  aromatic  proton  moieties  to 
be  determined  (Supporting  Information,  Table  S3).  Specifi¬ 
cally,  the  1H  signals  at  7.2  and  8.7  ppm  correspond  to  distinct 
thienothiophene  moieties  and  provide  important  quantitative 
constraints  that  allow  the  structure  of  the  bimolecular  crystal  to 
be  confidently  established.  A  detailed  discussion  of  these  anal¬ 
yses  are  beyond  the  scope  of  this  paper  and  will  be  published  in 
the  near  future. 

Considering  the  XRD  and  NMR  spectroscopy  results,  a  mon¬ 
oclinic  unit  cell  with  one  PBTTT-C14  monomer  per  fullerene 
molecule,  straight  side  chains,  and  planar  polymer  backbones 
is  a  plausible  starting  point  for  MM  modeling  of  the  structure. 
To  fine  tune  the  atomic  positions,  the  unit-cell  parameters  and 
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Figure  3.  A)  2D13C{1H}  HETCOR  NMR  spectrum  ofa  PBTTT-C1 6:PC71  BM 
biend  with  a  molar  ratio  of  1:1.  A  corresponding  1 D  13C{’H}  CP/MAS  spec¬ 
trum  is  shown  along  the  top  horizontal  axis,  and  a  1  D  single-pulse  1 H  MAS 
spectrum  is  shown  along  the  left  vertical  axis.  The  orange  and  blue  cir¬ 
cles  indicate  spectral  contributions  from  the  PBTTT-C16  and  the  PC71BM, 
respectively.  Lines  are  shown  to  aid  the  identification  of  correlated  signal 
intensities  that  establish  spatial  proximities  between  specific  moieties  of 
the  conjugated  polymer  and  the  fullerene  derivative.  The  red  arrows  mark 
correlated  signals  between  the  13C  moieties  of  the  PC71BM  and  the  'H  moie¬ 
ties  of  the  PBTTT-C16,  consistent  with  the  proposed  crystal  structure.  The 
light-blue  arrows  mark  correlated  signals  between  the  aromatic  13C  moieties 
of  the  conjugated  polymer  and  the  protons  from  the  PC71BM  side-chain, 
revealing  local  structural  disorder  of  fullerene  orientations  in  the  polymer 
crystal.  B)  Aliphatic  regions  of  1  D  13C{'  HjCP/MAS  spectra  of  the  PBTTT- 
C16:PC71BM  blend  (red)  and  the  neat  polymer  (black)  recorded  under  MAS 
conditions  of  12.5  kHz  with  a  CP  contact  time  of  2  ms.  The  colored  circles 
assign  the  13C  moieties  of  the  PBTTT-C16  side  chains  in  the  crystalline  (c) 
and  amorphous  (a)  regions  to  their  respective  signals  in  the  2D  spectrum. 

atomic  positions  were  relaxed  to  minimize  the  energy  associated 
with  intermolecular  interactions.  The  simulations  were  con¬ 
ducted  iteratively  until  there  was  good  agreement  of  the  measured 
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Figure  4.  An  extended  space-filling  structure  of  the  PBTTT-C1 4:PC71  BM 
bimolecular  crystal  based  on  the  simulated  unit  cell.  The  polymer  back¬ 
bones,  polymer  side  chains  and  fullerene  molecules  are  red,  orange  and 
blue,  respectively. 


XRD  patterns  and  the  local  proton  distances  determined  by 
1H  NMR  spectroscopy  (Supporting  Information,  Table  S3) 
with  those  calculated  for  the  minimum-energy  structure.  The 
resulting  structure,  which  is  shown  in  Figure  4  (and  is  provided 
in  cif  format  in  the  Supporting  Information),  has  twisted  and 
bent  polymer  backbones,  bent  side  chains,  and  is  triclinic  with 
the  angles  /3  and  yboth  equal  to  89°  (e.g.,  close  to  the  90°  ini¬ 
tial  guess).  A  1  ns  molecular  dynamics  (MD)  simulation  for  a 
2x3x2  super  cell  at  room  temperature  showed  stabilization 
both  in  energy  and  density  after  a  short  period  of  time  (200  ps). 
Both  the  MM  energy-minimized  structure  and  one  acquired 
from  a  snapshot  of  the  2x3x2  super  cell  obtained  towards 
the  end  of  the  MD  simulation  correlate  well  with  each  other  in 
terms  of  their  correspondence  with  the  proton-proton  distances 
determined  by  the  2D  NMR  spectroscopy  measurements  and 
in  their  use  to  simulate  the  2D  GIXS  pattern  (Supporting  Infor¬ 
mation,  Figure  S4). 

The  experimental  and  calculated  specular  XRD  intensities 
are  compared  in  Figure  2C.  These  reflection  intensities  are 
sensitive  to  the  electron  density  perpendicular  to  the  substrate. 
The  PC71BM  position  and  orientation,  the  polymer-backbone 
tilt  angle,  and  the  PC71BM  mole  fraction  primarily  determine 
the  electron  density  in  the  normal  direction.  The  experimental 
intensities  were  corrected  for  the  Lorentz  factor,  the  change  in 
irradiated  volume  with  q,  the  polarization,  and  the  Debye-Waller 
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factor  to  allow  comparison  with  intensities  calculated  from  the 
simulated  unit  cell.  The  good  agreement  between  the  experi¬ 
mental  intensities  (blue)  and  the  calculated  intensities  for  the 
simulated  unit  cell  (orange)  supports  the  accuracy  of  the  simu¬ 
lated  unit  cell.  In  particular,  the  model  captures  the  unexpect¬ 
edly  intense  (700)  reflection. 

Comparison  of  the  experimental  2D  GIXS  pattern  (Figure  2B) 
with  the  calculated  pattern  for  the  simulated  unit  cell  (Figure  2D) 
shows  good  agreement  of  the  reflection  positions.  Neverthe¬ 
less,  the  experimental  GIXS  pattern  exhibits  several  reflections 
that  are  broader,  compared  to  the  calculated  pattern,  which  are 
attributed  to  defects  (i.e.,  lamellar  stacking  disruptions)  and 
disorder  within  the  bimolecular  crystal  that  are  not  captured 
in  the  modeling.  The  MM  simulations  revealed  several  similar 
structures  with  comparable  total  energies.  Such  a  flat  energy 
landscape  allows  geometric  fluctuations  within  the  bimolecular 
crystal,  which  can  account  for  the  broadening  of  the  scattering 
reflections. 

The  solid-state  NMR  spectroscopy  results  also  confirm  dis¬ 
order  in  the  polymer  side  chains  and  the  fullerene  orientation. 
For  instance,  the  light  blue  arrows  in  the  2D  13C{1H}  FIETCOR 
spectrum  shown  in  Figure  3A  mark  the  correlated  intensity 
between  the  13C  signals  at  131,  135,  and  140  ppm  from  the 
aromatic  moieties  of  the  conjugated  polymer  backbone  and  the 
'H  signals  at  7.9  ppm  from  the  phenyl  moieties  of  the  PC71BM 
side  group.  Such  interactions  can  only  be  realized  with  a  dif¬ 
ferent  orientation  of  PC71BM  in  the  polymer  crystal,  reflecting 
a  distribution  of  local  configurations  that  are  not  apparent  in 
the  structure  shown  in  Figure  4  and  Figure  5.  In  this  case,  the 
PC71BM  moiety  must  be  flipped  by  180°,  such  that  the  PC71BM 
side  group  points  toward  the  PBTTT-C16  backbone.  Thus,  NMR 
spectroscopy  provides  direct  proof  of  local  structural  disorder 
regarding  the  fullerene  orientation,  in  line  with  the  observed 
flat  energy  landscape  in  MM  simulations  and  broadening  of  the 
XRD  reflections.  Such  orientational  disorder  may  cause  disrup¬ 
tions  of  the  ID  fullerene  channels  and  decrease  the  efficiency 
of  electron  transport  via  these  pathways. 

Additionally,  the  insertion  of  PC71BM  into  the  PBTTT- 
C16  crystal  results  in  local  disorder  of  the  alkyl  side  chains. 
Figure  3B  shows  the  aliphatic  regions  of  ID  13C{1H}  CP/ 
MAS  NMR  spectra  of  neat  PBTTT-C16  and  its  1:1  blend  with 
PC71BM.  Colored  circles  assign  the  signals  to  their  respective 
structural  moieties.  For  the  neat  polymer,  each  moiety  gives 
rise  to  two  spectral  contributions  with  slightly  different  13C  iso¬ 
tropic  chemical  shifts;  the  high-frequency  contributions  arise 
from  13C  moieties  with  y-neighbors  in  all-trans  conformations, 
a  necessary  prerequisite  for  crystalline  regions,  and  the  low-fre¬ 
quency  contributions  are  due  to  gauche  conformations,  which 
are  less  ordered  and  not  present  in  crystalline  domains. I22'23! 
For  neat  PBTTT-C16,  the  main  fraction  of  the  13C  signals  stems 
from  well-ordered,  crystalline  side  chains,  although  there  is 
some  intensity  from  disordered  side  chains,  which  manifests 
the  presence  of  a  significant  fraction  of  disordered  PBTTT- 
C16  side  chains.  Upon  intercalation  of  PC71BM,  almost  all  of 
the  aliphatic  13C  signal  intensity  in  the  PBTTT-C16:PC71BM 
blend  arises  from  such  disordered  side  chains.  Disorder  in 
the  PBTTT-C14  side  chains  in  the  bimolecular  crystal  was  also 
verified  by  observing  Fourier  transform  IR  (FTIR)  spectroscopy 
spectral  shifts  in  the  anti-symmetric  stretching  frequencies  of 


the  methylene  units  in  the  alkyl  chains  in  PBTTT-C14:PC71BM 
blends  as  a  function  of  composition  (Supporting  Information, 
Table  S5).  Thus,  the  good  agreement  between  the  calculated 
and  experimental  specular  XRD  intensities  and  2D  GIXS  pat¬ 
terns  suggests  that  the  structure  of  the  bimolecular  crystal 
closely  resembles  that  of  the  simulated  unit  cell,  while  the  2D 
NMR  spectroscopy  intensity  correlations  and  infrared  spectros¬ 
copy  show  some  variation  in  the  fullerene  orientation  and  side- 
chain  conformations  throughout  the  crystal. 

Expanded  structures  based  on  the  simulated  bimolecular- 
crystal  unit  cell  are  shown  in  Figure  4,5  with  the  polymer  back¬ 
bones,  polymer  side  chains,  and  fullerene  molecules  shown 
in  red,  orange,  and  blue,  respectively.  An  enlarged  version  of 
Figure  5D  with  a  stereoscopic  image  for  viewing  with  red/blue 
3D  glasses  is  shown  in  Figure  S9  in  the  Supporting  Information. 
These  figures  emphasize  important  features  of  the  bimolecular 
crystal  and  the  significant  differences  between  the  molecular 
packing  in  pure  PBTTT-C14I24,25)  and  the  PBTTT-C14:PC71BM 
bimolecular  crystal.  Such  features  are  similar  to  those  observed 
for  films  using  PBTTT  with  -C12H25  and  -C16H33  side  chains, 
which  exhibit  similar  XRD  patterns  (Supporting  Information, 
Figure  S5,  S6).  The  most-obvious  differences  between  the  struc¬ 
tures  are  the  increase  in  the  spacing  of  the  (100)  planes  due  to 
fullerene  insertion  (Figure  5A,D)  and  the  formation  of  fullerene 
channels  that  only  support  ID  electron  transport;  the  PBTTT- 
C14  side  chains  and  backbones  interfere  with  electron  trans¬ 
port  in  the  other  directions  (Figure  4,5E).  In  addition,  the  side 
chains,  which  tilt  at  an  angle  of  =35°  in  the  neat  PBTTT-C14 
polymer  (Figure  5A),  bend  in  the  bimolecular  crystal  to  allow 
dense  packing  in  the  interdigitated  regions,  while  maintaining 
enough  space  for  fullerene  intercalation  near  the  PBTTT-C14 
backbone  (Figure  5D).  Further  evidence  for  the  side-chain  con¬ 
formations  shown  in  Figure  5A,D  was  obtained  by  confirming 
that  the  lamellar  spacings  from  XRD  patterns  of  pure  PBTTT 
and  PBTTT:PC71BM  blends  with -C12FI25,  -Ci4H29,  and-C16FI33 
side  chains  match  those  predicted  by  the  modeling  (Supporting 
Information,  Figure  S5).  Likewise,  the  backbone  stagger  along 
the  c-axis  changes  to  accommodate  the  intercalated  fullerene  in 
the  bimolecular  crystal  (Figure  5C,F). 

Fullerene  intercalation  also  substantially  changes  the  stacking 
of  the  PBTTT-C14  backbones,  as  depicted  in  Figure  5B,E.  The 
positions  of  the  polymer  backbones  alternate  up  and  down,  as 
half  of  the  backbones  interact  with  the  fullerenes  in  the  upper 
part  of  the  unit  cell,  while  the  other  half  interact  with  fuller¬ 
enes  in  the  lower  part  of  the  unit  cell  (Figure  5E).  Moreover,  the 
thienothiophene  units  closest  to  the  fullerene  molecules  twist 
with  respect  to  the  rest  of  the  polymer  backbone,  presumably 
to  improve  the  intermolecular  interactions  between  the  PBTTT- 
C14  backbone  and  the  fullerenes.  The  3D  renderings  of  the 
structures  (Supporting  Information,  Figure  S9)  show  that  the 
backbones  also  bend  slightly  and  wrap  around  the  fullerenes. 
This  twisting  and  bending  results  in  a  non-planar  polymer 
backbone  with  twist  angles  between  the  thienothiophene  and 
thiophene  moieties  ranging  from  10°  to  40°. 

These  molecular  packing  features  explain  important  prop¬ 
erties  of  PBTTT-C14:PC71BM  blends.  For  instance,  the  back¬ 
bone  twists  and  backbone-fullerene  interactions  explain  why 
the  PBTTT-C14:PC71BM  absorption  spectrum  differs  signifi¬ 
cantly  from  the  superposition  of  the  PBTTT-C14  and  PC71BM 
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Figure  5.  Molecular  structures  of  the  pure  PBTTT-C14  (A-C)  and  the  PBTTT-C14:PC71BM  (D-E)  bimolecular  crystal  viewed  along  the  b-axis  (A,D)  the 
c-axis  (B,E),  and  the  side  chains  (C,F).  The  polymer  backbones,  polymer  side  chains  and  fullerenes  are  red,  orange  and  blue,  respectively.  The  projection 
of  the  b  and  c  axes  and  the  angle  aare  shown  in  (C)  and  (F).  Note  that  these  structures  do  not  capture  the  disorder  present  in  the  films,  since  every  unit 
cell  in  these  structures  is  the  same.  Warren-Averbach  analysis  ofthe  specular  diffraction  peaks  gives  an  out-of-plane  crystallite  size  and  paracrystallinity 
parameter  of  =24  nm  and  =2%,  respectively.  I45,46l  Thus,  the  actual  structure  ofthe  films  will  contain  more  disorder  than  is  apparent  in  these  figures. 


absorption  spectra.!17,26'  Likewise,  the  previously  observed 
order-of-magnitude  decrease  in  the  PBTTT-C14  hole  mobility 
upon  the  addition  of  PC71BM  can  be  attributed  not  only  to  the 
dilution  of  the  hole-transporting  material,  but  also  to  the  dis¬ 
ruption  of  the  polymer  7l-stacking  that  results  both  from  the 
backbone  twists  and  the  additional  unit-cell  volume  required 
to  account  for  the  intercalated  fullerene.!10,27!  Indeed,  density 
functional  theory  (DFT)  calculations  (Supporting  Information, 
Table  S4)  reveal  that  the  electronic  coupling  between  polymer 
chains  along  the  backbone-stacking  direction  (fa-axis)  drops  by 
an  order-of-magnitude  in  the  presence  of  intercalated  PC7iBM. 


In  addition,  DFT  calculations  indicate  that  there  is  non-neg- 
ligible  electronic  coupling  (=20-40  meV)  for  electrons  between 
neighboring  PC71BM  molecules  along  the  fa-axis  in  the  bimo¬ 
lecular  crystal,  but  negligible  electronic  coupling  strength 
(=0  meV)  between  PC71BM  along  other  directions.  This  implies 
a  possible  electron-transport  pathway  through  the  ID  PC71BM 
channels.  The  very  low  electron  mobility  in  bimolecular-crystal 
films  can  be  attributed  to  the  ID  character  of  the  fullerene 
channels  and  the  disorder  of  the  fullerenes  in  the  bimolecular 
crystal.!10!  In  a  polycrystalline  film,  it  is  expected  to  be  difficult 
for  electrons  to  move  between  differently  oriented  crystals, 
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because  the  ID  electron  channels  will  not  generally  be  aligned. 
Moreover,  electron  transport  perpendicular  to  the  substrate, 
which  is  critical  for  solar  cell  operation,  will  be  especially  diffi¬ 
cult,  because  the  PBTTT-C14:PC71BM  bimolecular  crystals  gen¬ 
erally  form  lamellar  structures  with  the  electron  channels  run¬ 
ning  parallel,  rather  than  perpendicular,  to  the  substrate.  The 
low  electron  mobility  does  not  appear  to  be  a  result  of  a  large 
fullerene-fullerene  distance  in  the  bimolecular  crystal,  because 
the  average  fullerene  separation  in  the  simulated  unit  cell  is 
9.8  A,  which  is  comparable  to  the  reported  center-to-center  dis¬ 
tance  of  9.88  A  in  a  C70  crystal.'28'  On-going  quantum-chemical 
calculations,  based  on  the  optimized  structures,  should  prove 
useful  to  understand  the  electronic  properties  of  the  bimo¬ 
lecular  crystals  at  a  deeper  level. 

Since  fullerenes  intercalate  between  the  side  chains  of  many 
polymers  used  in  solar  cells,  PBTTT:PC71BM  is  a  valuable 
model  system.  While  it  will  not  be  possible  to  use  XRD  to  deter¬ 
mine  the  structure  of  other  polymer-fullerene  mixed  phases 
in  cases  where  there  are  only  a  few  XRD  peaks,  it  should  be 
possible  to  generate  approximate  structures  by  assuming  that 
fullerenes  alter  the  polymer  packing  in  ways  that  are  similar 
to  those  reported  in  this  study.  Molecular  mechanics  and 
molecular  dynamics  simulations  can  be  used  to  fine  tune  these 
structures  to  be  consistent  with  constraints  provided  by  NMR 
and  FTIR  spectroscopy.  Now  that  the  structure  of  the  impor¬ 
tant  PBTTT-C14:PC71BM  model  system  is  known,  it  will  be 
fruitful  to  compare  theoretical  property  predictions  based  on 
the  structure  with  measurements  made  by  powerful  spectro¬ 
scopic  techniques,  such  as  UV-vis,  NMR,  and  Raman  spec¬ 
troscopy.  When  the  accuracies  and  agreement  of  the  theoret¬ 
ical  and  spectroscopy  tools  are  validated,  it  will  be  possible  to 
use  them  with  less  ordered  materials  to  gain  information  that 
has  not  been  previously  obtainable.  We  expect  that  solid-state 
NMR  spectroscopy,  which  has  been  rarely  used  to  investigate 
bulk  heterojunctions, I21-29-31!  will  be  provide  helpful  insights  for 
determining  how  fullerenes  interact  with  the  complex  donor- 
acceptor  polymers  that  enable  the  highest  power  conversion 
efficiencies.  Knowing  more  about  how  molecules  pack  in  bulk 
heterojunctions  will  be  useful  to  researchers  trying  to  under¬ 
stand  why  recombination  is  much  slower  in  some  polymer- 
fullerene  blends  than  others  and  will  enable  the  design  of  more 
efficient  solar  cells.  More  generally,  the  complicated  interplay 
of  molecular  order  and  disorder  in  the  polymer:fullerene  bimo¬ 
lecular  crystal  is  intrinsic  to  numerous  other  semicrystalline 
materials,  especially  those  for  which  molecular  interactions  at 
surfaces  broadly  govern  macroscopic  properties.  These  include 
diverse  semicrystalline  polymers  or  hybrid  nanocomposites  that 
are  used  broadly  as  ion-conducting  membranes  in  fuel  cells  or 
batteries,'32'  as  environmental  barriers  or  for  structural  pack¬ 
aging  in  microelectronics,!33'  or  as  dye  hosts  in  opto-electronic 
applications.'34!  For  these  and  related  systems,  feedback  among 
complementary  synthesis,  characterization,  and  modeling  strat¬ 
egies  are  key  to  improving  material  properties  to  reach  device 
and  technology  performance  targets. 

Experimental  Section 

Materials:  PBTTT  was  synthesized  with  -C12H25,  -C14H29  and  -Ci6H33 
side  chains  as  described  by  McCulloch  et  all27!  The  number-average 


molecular  masses  (M„)  of  PBTTT  with  -G|2H25,  -C3  4H29  and 
-C16H33  side  chains  were  determined  to  be  28,  22,  and  35  kDa  with 
polydispersities  of  1.8,  2.0,  and  2.0,  respectively,  using  gel-permeation 
chromatography  (GPC)  against  polystyrene  standards.  PC71BM  was 
purchased  from  NanoC.  The  PBTTT:PC71BM  films  for  X-ray  diffraction 
were  spin-cast  from  ortho-dichlorobenzene  on  octadecyltrichlorosilane 
(OTS)-coated  silicon  substrates  and  annealed  at  180  °C  for  10  min. 
For  the  NMR  spectroscopy  measurements,  PBTTT-C1 6:PC71  BM  films 
with  a  molar  ratio  of  1:1  were  drop-cast  from  ortho- dichlorobenzene  on 
octadecyltrichlorosilane  (OTS)-coated  silicon  substrates  and  annealed  at 
180  °C  for  10  min.  The  structural  similarity  of  drop-cast  and  spin-cast 
films  was  confirmed  by  XRD. 

X-Ray  Diffraction:  2D  GIXS  measurements  were  performed  at  the 
Stanford  Synchrotron  Radiation  Laboratory  (SSRL)  on  beamline  11-3  at 
an  energy  of  12.7  keV  and  an  incidence  angle  of  0.10°  using  a  MAR345 
image  plate  area  detector.  We  omitted  the  features  near  =  0  from  our 
experimental  2D  GIXS  data,  as  these  peaks  are  not  accurately  measured 
in  image  plate  collection  geometry.  I35!  High-resolution  specular  diffraction 
and  high-resolution  grazing-incidence  X-ray  diffraction  were  performed 
at  an  energy  of  8  keV  on  SSRL  beamlines  2-1  and  7-2,  respectively.  To 
produce  PBTTT-C14:PC71BM  films  with  a  preferred  in-plane  orientation, 
we  aligned  pure  PBTTT-C14  using  a  flow-coating  method'18'  and 
subsequently  spin-cast  a  PC71BM  layer  on  top  of  the  aligned  PBTTT-C14 
film  by  using  the  orthogonal  solvent  dichloromethane  (DCM).  Annealing 
allowed  diffusion  between  the  layers  and  resulted  in  the  coexistence  of 
aligned  PBTTT-C14,  aligned  bimolecular  crystals  and  randomly-oriented 
PC71BM  as  shown  in  the  specular  X-ray  diffraction  data  (Supporting 
Information,  Figure  SI). 

IR  Spectroscopy:  IR  spectroscopy  was  performed  using  a 
commercial  FTIR  spectrometer.  The  samples  were  mounted  on  a 
custom  transmission  stage  such  that  the  IR  beam  was  incident  on  the 
samples  nominally  at  Brewster’s  angle  (about  72°).  The  polarization 
of  the  IR  beam  was  controlled  by  a  wire-grid  polarizer  placed 
immediately  before  the  sample  stage.  The  IR  spectroscopy  results 
were  fitted  with  a  least-squares  regression  of  both  s-  and  p-polarized 
spectra  simultaneously. 

Theoretical  Modeling:  Molecular  mechanics  (MM)  and  molecular 
dynamics  (MD)  simulations  were  performed  on  the  PBTTT-C1 4:PC71  BM 
unit  cell  using  the  universal  force  field,  as  described  by  Rappe  et  al.,l36l 
as  employed  in  the  Materials  Studio  software  suite.  The  SimDiffraction 
code'37'  was  used  to  simulate  the  XRD  patterns  of  the  optimized  blend 
unit  cells.  The  a*  direction  of  the  crystal  unit  cell  was  set  to  be  oriented 
parallel  to  the  substrate  normal.  The  diffraction  simulation  conditions 
were  set  to  be  identical  to  the  experiment.  The  in-plane  and  out-of- 
plane  orientations  were  set  as  infinity  and  10°,  respectively,  to  describe 
the  preferential  orientation  distribution  of  the  films.  The  intensity  of  the 
diffracted  peaks  was  calculated  by  accounting  for  the  Lorentz  factor, 
the  polarization  factor,  and  the  Debye-Waller  factor,  as  described  in  the 
literature.'37!  Density  functional  theory  calculations  were  carried  out 
at  the  B3LYP/6-31G**  level.  The  evaluation  of  the  nearest-neighbor 
effective  transfer  integral  (electronic  coupling)  was  determined  through 
the  use  of  a  fragment  orbital  approach  in  combination  with  a  basis  set 
orthogonalization  procedure.'38' 

Solid-State  NMR  Spectroscopy:  Solid-state  NMR  spectra  were 
recorded  on  a  Bruker  AVANCE  II  NMR  spectrometer  with  an  11.7  Tesla 
wide-bore  superconducting  magnet,  operating  at  Larmor  frequencies 
of  500.24  MHz  for  'H  and  125.79  MHz  for  13C.  The  experiments  were 
conducted  at  293  l<  on  films  that  were  scraped  off  substrates  and  loaded 
into  4.0  mm  zirconia  rotors  with  Kel-F  caps  and  measured  using  a 
4.0  mm  double-resonance  H-X  magic-angle-spinning  (MAS)  probehead 
under  MAS  conditions  at  a  rotation  rate  of  12.5  kHz.  The  spectra 
were  acquired  with  a  'H  90°-pulse  length  of  2.5  |j.s,  corresponding  to 
a  radio  frequency  (rf)-field  strength  of  100  kHz,  and  a  relaxation  delay 
of  1.5  s.  For  the  ID  13C  cross-polarization  (CP)/MAS  and  2D  13C{'H} 
HETCOR  spectra,  magnetization  was  transferred  from  'H  to  13C  nuclei 
through  CP  by  adiabatic  passage  with  a  contact  time  of  2  ms. I39!  A 
broadband  SPINAL-64  heteronuclear  'H  decoupling  scheme  was  used 
during  detection.!40'  For  the  2D  ^CfH}  HETCOR  spectrum,  quadrature 
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detection  in  the  indirect  'H  dimension  was  achieved  by  using  the 
States-TPPI  detection  scheme.!41!  Homonuclear  dipolar  couplings 
affecting  the  'H  dimension  were  reduced  by  applying  high  power  'H- 
'H  homonuclear  decoupling  using  the  eDUMBO-l22  sequence  during 
the  proton  spin  evolution  time  t  j  with  a  phase-modulated  rf  pulse  of 
constant  amplitude  (100  kHz). I42!  80  t,  increments  of  128  ps  with  848 
transients  per  increments  were  recorded.  In  the  indirect  'H  dimension, 
a  scaling  factor  of  0.601  was  determined  experimentally  with  a 
2D  1 2 3 H {4  H}  homonuclear  correlation  experiment.  The  chemical  shifts  are 
reported  relative  to  tetramethylsilane  (TMS)  using  tetrakis (tri methylsi lyl) 
silane  (TKS)  and  adamantane  as  secondary  standards  for  'H  and  13C, 
respectively.!4344! 
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